Confinement induced resonances in anharmonic waveguides 
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We develop the theory of anharmonic confinement-induced resonances (ACIR). These are caused 
by anharmonic excitation of the transverse motion of the center of mass (COM) of two bound 
atoms in a waveguide. As the transverse confinement becomes anisotropic, we find that the COM 
resonant solutions split for a quasi-lD system, in agreement with recent experiments. This is 
not found in harmonic confinement theories. A new resonance appears for repulsive couplings 
(asB > 0) for a quasi-2D system, which is also not seen with harmonic confinement. After inclusion 
of anharmonic energy corrections within perturbation theory, we find that these ACIR resonances 
agree extremely well with anomalous ID and 2D confinement induced resonance positions observed 
in recent experiments. Multiple even and odd order transverse ACIR resonances are identified in 
experimental data, including up to A'^ = 4 transverse COM quantum numbers. 
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I. INTRODUCTION 



Ultra-cold low-dimensional atomic gases show unique 
quantum properties, and have attracted a great deal 
of interest. For a one-dimensional (ID) Bose gasfl- 
|3|, the finite-temperature correlations predicted for a 
To nks- Girardeau gas0, Q have been experimentally 
verified[6, 7], and a cross-over to a non-equilibrium super 
To nks- Girardeau gas has been realized [8,]. For a two- 
dimensional (2D) geometry, the Berezinskii-Kosterlitz- 
Thouless phase transition was predicted P, \1^ , and sub- 
sequently observed in experiment [Il| . In these experi- 
ments, one or two spatial degrees of freedom are removed 
by introducing tight confinement via an optical lattice or 
a tightly focused anisotropic dipole trap. 

Atomic interactions can also be tuned precisely by 
means of a molecular Feshbach resonance in an external 
magnetic field (12h14| . This allows an effective contact 
interaction with scattering length a^o to be created, with 
scattering lengths that can be varied over a wide range 
of positive and negative values. Owing to these methods, 
low-dimensional atomic gases now provide a high degree 
of control for tests of fundamental many-body physics in 
reduced dimensions. These systems are often much sim- 
pler than condensed matter physics experiments, which 
have complex crystal structure, interactions and disorder. 

Confinement-induced resonances (CfR) are one of the 
most intriguing phenomena found in low-dimensional sys- 
tems. These were first predicted theoretically by 01- 
shanii |l3| , who considered a two-body S-wave scattering 
problem in a quasi- ID trap with cylindrically symmetric 
transverse harmonic confinement. The CIR can be un- 
derstood as a novel type of Feshbach resonance, where 
the transverse ground mode and the manifold of molec- 
ular internally excited modes play the roles of the open 
and closed channels respectively \l&. Related effects oc- 
cur in mixed dimensional traps |l7l . Il8| . A direct gen- 
eralization of Olshanii's theory to anisotropic transverse 
confinement shows that there is only one harmonic CIR 
(HCfR) resonance, no matter how large the transverse 



anisotropy [I^. For large anisotropy, this theory crosses 
over smoothly to the case of a quasi-2D trap, where a 
single HCIR occurs with a negative S-wave scattering 
length, asD < ^^^. 

There have been a number of related experimental in- 
vestigations, which in some cases appear to contradict 
each other. In the recent Innsbruck Cs experiment with 
a quasi- ID geometry [23| . two or more resonances were 
observed as the transverse confinement became more and 
more anisotropic. For a quasi-2D geometry, some exper- 
iments have observed 2D resonances on the attractive 
side with a^o < [2J|, while others have resonances on 
the repulsive side with a^o > |23l . [25j . Both the ob- 
servation of multiple resonances and 2D resonances with 
repulsive interactions are in disagreement with standard 
HCIR predictions ^^. 

In this paper, a detailed explanation is proposed for 
these anomalous resonances. The mechanism is that the 
new resonances are due to center-of-mass (COM) exci- 
tations of molecules or atom pairs. These have a differ- 
ent character to the excitation of internal molecular de- 
grees of freedom found in the Olshanii approach and its 
generalizations. The new COM resonances can only be- 
come coupled to the input state by anharmonic terms in 
the trapping potential. Hence, we term these effects an- 
harmonic confinement induced resonances (ACIR). The 
ACIR resonances cannot occur in harmonic traps due to 
Kohn's theorem [26]. However, they are certainly observ- 
able in current ultra-cold atomic physics experiments, 
which have relatively large anharmonicities. 

The coupling of the COM motion to the relative mo- 
tion gives additional degrees of freedom not found with 
parabolic traps. This causes a series of additional scat- 
tering resonances due to the mixing of COM and rel- 
ative motion. The nonlinear mixing caused by anhar- 
monic terms in the potential makes these phenomena 
analogous to frequency-mixing effects found in nonlinear 
optics. They provide a fundamentally new pairing mech- 
anism, which may lead to new opportunities for quantum 
engineering in atomic, photonic or acoustic waveguides. 
Our results are therefore qualitatively different to har- 



monic CIR. We predict both multiple ID resonances and 
2D resonances with repulsive interactions. Both results 
are in quantitative agreement with experiment. 

We note that this possibility was also envisaged in 
three earlier papers. In Peano et al |27i], the general idea 
is addressed, although using a different technique, and 
for a different type of trap. Kestner and Duan [28] treat 
anharmonic resonance in a double well. In a more recent 
investigation, parallel to our own, Sala et al [2^ have also 
concluded that the recent Innsbruck experiments provide 
evidence for ACIR resonances. The main differences in 
the treatment are that we have accurately calculated the 
size of the anharmonic resonance shifts, as well as giv- 
ing quantitative estimates of relevant parameters. We 
also compare our theory with the observed multiple res- 
onances, including even and odd order COM resonances. 

The paper is arranged as follows. We first analyze the 
types of transverse excitations available, and the opera- 
tional processes that can lead to the observed resonances 
(Sec|lT|. In Sec lIIIl a Hamiltonian model of nonlinear 
CIR is presented, by introducing anharmonic perturba- 
tions in the Hamiltonian. In Sec lIVl this is analyzed us- 
ing perturbation theory for the ID case. The results are 
compared to experiments on anisotropic traps, showing 
the observed splitting is well- explained with the anhar- 
monic COM resonance ACIR model. Next, we consider 
results for the case of large trap anisotropics, and demon- 
strate that the observed resonances can be quantitatively 
explained with excellent accuracy by considering multi- 
ple resonances with both even and odd order COM trans- 
verse quantum numbers. A similar calculation is carried 
out for a quasi-2D system in Sec. |Vl which is also com- 
pared with experiment. The main results are summarized 
in SecED 



II. TWO-BODY CIR PHYSICS 

In recent one and two dimensional confinement induced 
resonance experiments, there are many observed reso- 
nances not explained by conventional CIR theory. In 
two dimensions, resonances are observed for a^o > 0, 
(the BEC side of the resonance) which is the opposite to 
that expected in the usual theory. Similarly, unexplained 
multiple resonances occur for one-dimensional CIR with 
anisotropic transverse confinement. These are also not 
predicted by the simple two-particle model [15] with lin- 
ear confinement. However, the experiments have some 
features not included in this idealized model, and the 
obvious question is: which experimental properties are 
responsible for the additional observed resonances? 



A. Harmonic CIR solutions 

The possible modes of excitation of a pair of atoms in 
a transverse potential are illustrated in Fig[TJ The COM 
mode is shown in (a), with two atoms moving together. 





Figure 1. (Color online) (a) Illustration of an ACIR reso- 
nance, in which two atoms are driven in a molecule-like ex- 
cited state with a transverse COM motion in the trap. Both 
atoms move together, to give a transverse resonance of the 
COM degree of freedom. This is only coupled to the atomic 
ground state if the potential is anharmonic. (b) Illustration 
of an HCIR resonance, in which the atoms move in relative 
motion, to give an excited state rather like a conventional 
vibrational excited state of a molecule. 



This is not coupled to the atomic ground state in a har- 
monic trap, due to Kohn's theorem. The relative motion 
mode is shown in (b) , with an excitation of the relative 
coordinate. This is the usual harmonic confinement in- 
duced resonance in a one dimensional waveguide with a 
transverse parabolic potential. The HCIR resonance is 
simply the first internally excited resonance of the two- 
body ground state. 

However, there is a subtlety here. Both types of excita- 
tion are an adiabatic continuation of the transversely ex- 
cited free-particle states, as the inter-particle interaction 
is increased. The free-particle states have both center- 
of-mass (COM) quantum numbers N^^Ny, and internal 
quantum numbers, nx,ny. Therefore, there are a large 
number of possible excited states that could, in principle, 
be coupled to the atomic ground state in a Feshbach type 
of resonance. 

In a rotationally symmetric model, both the internal 
and COM excited states near 2hiij are degenerate in en- 
ergy, and cannot lead to distinct resonances even if cou- 
pled to the incoming states. If there is no rotational 
symmetry, this degeneracy is broken, leading to the pos- 
sibility of multiple resonances. Experimentally there is 
more than one resonance in the recent Innsbruck exper- 



imcnt f2Jy|. Due to the similarity with the HCIR predic- 
tions [15J , these were initiaUy identified as conventional 
harmonic CIR resonances, corresponding to excitations 
of internal degrees of freedom. 

However, an analysis of the anisotropic case [l^ indi- 
cates that internally excited two-body states only give 
rise to only a single CIR. We note that Kohn's theorem 
|2a | prohibits the excitation of COM resonances from in- 
coming states with no transverse momentum, if the con- 
finement is parabolic. One may ask why there are not 
multiple CIR eigenstates with different internal energies 
due to the different confinement strength in orthogonal 
trap directions. The reason for this is due to the singular 
nature of the interaction. 

Consider what happens to a single particle in a rota- 
tionally symmetric potential, which corresponds to the 
internal quantum numbers in a relative coordinate pic- 
ture of a two-particle problem. In a 2D system of nonin- 
teracting particles in a harmonic oscillator potential, the 
ground state has n^ — ny ^ 0. The internally excited 
states can therefore be labelled either by their internal 
harmonic oscillator quantum numbers for relative mo- 
tion, or by their angular momentum quantum numbers. 
This, one can have Ux = riy = 1, or else a radial quantum 
number rir and a magnetic quantum number m. Due to 
the singular potential at the origin, only S-wave incoming 
states with m — experience any coupling. This leads to 
a relatively low-lying excited state due to the coupling, 
and hence to a single CIR. This state is adiabatically de- 
formed when the symmetry is broken, without leading to 
a second CIR. 

This single degenerate CIR exists on the positive side 
of the Feshbach resonance i.e, a^D > for a quasi-lD 
system. It transfers to the negative side, i.e., a^n < 
for a large asymmetry or a quasi-2D system. This last 
conclusion is compatible with other calculations of 2D 
CIR. However, these conclusions only take into account 
the internal energy of a two-particle state, not the COM 
energy. 



there are more degrees of freedom available for excitation, 
since the COM quantum numbers must now be included 
in the description. 

The interesting issue is whether these observed CIR 
effects can be explained as anharmonic (ACIR) res- 
onances due to center-of-mass excitations of resonant 
bound states [27l-(2l|. This is iUustrated in FigdJ by the 
two atoms moving together in (a). Such effects can only 
occur in an anharmonic trapping environment, which al- 
lows coupling between incoming scattering states with 
zero transverse excitation and an outgoing transverse 
COM excitation. 

In another 2D experiment on ultracold ^°X at 
Cambridge|24j. a CIR resonance occurs on the attractive 
side of a Feshbach resonance, as expected. This exper- 
iment has a much lower anharmonicity than the Inns- 
bruck experiment. It also uses a different technique to 
identify the resonance, employing RF spectroscopy rather 
than molecular losses. Thus, there are distinct resonance 
signatures used in the two published experiments. The 
Cambridge data appears to show evidence for anomalous 
resonance features on the a^D > side of the Feshbach 
resonance, but this effect is greatly reduced compared to 
the Innsbruck observations. 

All Feshbach bound states have a bound molecular 
fraction[30] in which the atoms have a small separation. 
We conjecture that this molecular fraction is larger for 
COM ACIR resonances (a), where the atoms are in a rel- 
ative ground state - compared to internally excited HCIR 
resonances (b), where the atoms are in a relative excited 
state. This would mean that COM excitations would 
have a relatively larger 3-body recombination loss due 
to molecule formation. This is precisely the signature of 
the resonances used in the Innsbruck experiments. RF 
spectroscopy, used in the Cambridge experiments, on the 
other hand has different characteristics. Thus, it is not 
unreasonable to expect the two experiments to have a dif- 
ferent relative sensitivity to internal and COM molecular 



resonances. 



B. Anomalous CIR experiments 

In recent bosonic experiments on ultracold ^^"^Cs at 
Innsbruck[23| , a strong, transversely anisotropic quasi- 
lD confinement is used to create an initially strongly re- 
pulsive {asD > 0) Tonks gas. This is followed by a sudden 
change in i?-field to a new value, resulting in a molecu- 
lar loss signature for a resonance which is confinement 
dependent. Numerous multiple confinement-induced res- 
onances are observed. There is even an unexpected res- 
onance for a^D > in the two-dimensional (2D) limit, 
which has also been measured using release energy data 
in fermionic ^Li experiments [25|. All these observations 
contradict the harmonic waveguide theory given above. 

However, it is important to recognize that the waveg- 
uide potential in these experiments are generally anhar- 
monic, so that Kohn's theorem does not apply. Hence, 



III. ANHARMONIC WAVEGUIDES 

For technical reasons explained below, current CIR ex- 
periments typically involve an anharmonic confinement 
mechanism. In such cases, the excitation of a COM de- 
gree of freedom can be coupled to input states with no 
transverse COM excitation. This coupling mechanism 
would explain observed anomalies such as multiple res- 
onances, different to those predicted using the standard 
parabolic confinement theory. 

There are three possibilities which might allow reso- 
nant coupling to additional confinement induced bound- 
states through nonlinear mechanisms: 

1. The Kohn theorem only applies for parabolic con- 
finement. Experimentally the optical confinement 
is sinusoidal and/or Gaussian, not parabolic. This 
allows direct coupling to N^ = 2 states, or even 



Nsc = 1 states, depending on the type of anhar- 
monicity. 

2. In some experiments the width of the Feshbach 
resonance is comparable to the separation of the 
transverse modes, ahowing contributions from the 
molecular bound state channel as well as the atomic 
channel. This still requires anharmonic coupling to 
access bound states having COM transverse energy. 

3. At high density the mean- field background poten- 
tials of the other atoms may provide an anharmonic 
effective potential which is not parabolic. In such 
cases, one may expect the collective oscillation fre- 
quencies to play a role. 

Coupling to COM excitations is possible whenever the 
transverse response is nonlinear, and the potential de- 
parts from a parabolic shape. This is not inconsistent 
with the ultra-cold atomic physics experiments, which 
generally involve sinusoidal laser trapping potentials. 
These are only approximately parabolic in the strongly 
confined limit. It is this possibility of nonlinear CIR due 
to anharmonic confinement which is explained below. Al- 
though related theoretical work that has been carried out 
includes one or the other of these effects, it is important 
to include both anisotropy and anharmonicity to fully ex- 
plain the observed ACIR resonances. 



A. Hamiltonian 

In this paper, we consider the simplest model of non- 
linear CIR, with a single-channel S-wave interatomic po- 
tential and an anharmonic trapping potential. We do 
not take into account either many-body corrections or 
explicit molecule formation channels [30]. This model is 
therefore applicable to relatively dilute quantum gases 
with a broad Feshbach resonance. To model the non- 
linear CIR effect in greater detail, consider two atoms 
with mass m which are anisotropically confined in the 
transverse direction, and can almost freely move in the z 
direction. Such a model can treat both a quasi-lD and 
quasi-2D experiment, by taking one of the trapping fre- 
quencies to zero. The Hamiltonian of the two atoms in a 
quasi-lD system is, therefore: 



where 



H = Hi+H2 + U{ri-r2) , 



U{r) =53c(5(r)— r 



(3.1) 



(3.2) 



is the atomic interaction for S-wave scattering de- 
scribed by a zero-range pseudo-potential with interaction 
strength g^jj = ATrh^a^D/m corresponding to a scatter- 
ing length of a3£), and Hi is a single-particle Hamiltonian 
including external potential and kinetic energy terms. 



1. Anharmonic confinement potential 

For dipole trapping experiments with ID and 2D opti- 
cal lattices, the trapping potential near a potential min- 
imum at r = is due to an optical standing wave of 
form: 
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(3.3) 



Here, Vx.y (r) are the two orthogonal slowly varying 
potential energy envelopes of standing waves due to the 
atomic dipole interactions with the two trapping lasers 
at optical wavelengths A^;, Aj,. Thus, V^° are potential 
well-depths, leading to trap frequencies tu^, tUy in the x, y 
directions. We have used a scale length of the reduced 
oscillator lengths dx^y = y''2h/mujx.y in each direction. 
It is also common to use the single atom oscillator length 
definition of ax,y = y^h/mujx,y, which we will utilize in 
comparisons with experiment in later sections. 

To next order beyond the linear confinement approxi- 



mation, we have introduced ax,cty 
anharmonic parameters, so that: 



<C 1 as the dominant 
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(3.4) 



For plane waves, these quartic anharmonic terms are 
the lowest order possible. More generally, the potential 
may be neither parabolic nor sinusoidal. Examples of 
this include the potential found in an optical fibre, which 
can be engineered to any desired shape, and potentials 
found in experiments using magnetic trapping or focused 
Gaussian beams. For this reason, we expect cubic, quar- 
tic and higher order anharmonic parameters in any real 
experiment. However, the quartic term given above is 
due to spatial modulation on optical wavelength scales. 
This is generally larger than cubic anharmonic terms like 
r • WVx caused by focusing effects. 

For simplicity, we suppose that the the dominant an- 
harmonic effects are caused by anharmonic parameters 
ax,y, and the single-particle Hamiltonian is: 



H, 



2m 
1 

2 



1 



{l + axx^/dl) + 
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{I + ayyf /dl) .{i= 1,2) (3.5) 



Thus, the trapping Hamiltonian p. II) has the form of a 
harmonic Hamiltonian Hh plus an anharmonic terms iJf 
and Ry in the xernd y directions respectively: 

Hi+H2 = Hh + H: + HI 

= Hh + Ha . (3.6) 



Experiment 


^^^Cs 


40^ 


Trap frequency u] 


27r X U.5kHz 


27r X 80 kHz 


Wavelength A 


1.064 X 10^'^m 


1.064 X lO^^'m 


Atomic mass m 


2.22 X IQ-'^^kg 


0.6635 X 10"^^fcg 


Length dx,y 


0.102 X 10"'^m 


0.08 X 10"^m 


Anharmonicity ax,y 


-0.121 


-0.075 



Table I. Typical anharmonic parameters for CIR experiments 
using optical lattices, following data from Innsbruck (Cs) [23 1 
and Cambridge (K) [2J|. Quantitative values depend on the 
trap frequency, which is varied over a range of values. 



Next, we can estimate typical parameter values in re- 
cent experiments, as shown in Table ID We note that 
in the Innsbruck experiments with relatively large ob- 
served anomalies, the dimensionless anharmonic param- 
eter was typically 12%, which is substantially larger than 
in the case of the Cambridge experiment, with an anhar- 
monicity of 7.5%. These parameters are calculated for 
alkali metal atoms, micron wavelength lasers and typical 
10 — lOOkHz trap frequencies. Obviously, large changes 
in anharmonicities are easily obtained by changing any 
of the relevant factors. 

These anharmonic parameters lead to energy shifts in 
the ACIR resonances, which we calculate below. More 
importantly, any type of anharmonic potential allows a 
coupling between relative and COM motion, which is oth- 
erwise prohibited due to the Kohn theorem. 



B. COM energies 

We can now make a preliminary estimate of atomic and 
molecular energies in the two-particle sector, given the 
anharmonic trapping potential. These estimates assume 
sufBciently tight internal binding so that only the COM 
energies are changed by the anharmonicity. While this 
is not accurate near threshold, it allows an estimate of 
the size of anharmonic perturbation energies. It will also 
be used to check the validity of subsequent results in the 
tight-binding limit. 



1. Atomic energy 

For a single particle, the perturbation theory solution 
including the anharmonic parameter is well known. We 
can calculate how the free atomic energy. 
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is changed by the anharmonic perturbation. To first 
order in perturbation theory, the modified transverse 
ground state energy is an elementary perturbation theory 



result, such that Ea — E)^ ' + E)^ ' , where 
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Here, lipn^.ny) is the single-particle eigenstate of H^ , 
which is treated as the zero-order wave-function. For a 
threshold resonance experiment, the incoming total en- 
ergy of two atoms initially in a transverse and longitudi- 
nal ground state is therefore: 
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Using the numbers in Table (1), this indicates that 
the scale of anharmonic energy perturbations should be 
around 5 — 10% of the transverse trap frequency, for the 
parameters of recent experiments. 



2. Molecular energy 

In the COM relative-coordinate frame, with R = 
(ri -)- r2) /2 , r = ri — r2 , the harmonic term is: 



Hh — H, 
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and the anharmonic term is: 
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Here, r = (a;,y, z) , R = {X,Y,Z) , and M = 2m , 
/i ~ to/2 are the mass of the COM and the reduced 
mass, respectively. 

As a point of reference for the more detailed calcula- 
tions given below, we now consider how the COM energy 
of molecular bound states changes due to anharmonicity. 
For a broad Feshbach resonance, with strong binding so 
that the internal molecular energy is not changed by the 
waveguide, the internal molecular bound state energy in 
free space for an attractive interaction is known to be: 
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3D 



For a tightly bound molecule described only by the COM 
coordinates {X,Y), the oscillator frequencies are cux, ujy 
as before. To first order in perturbation theory, the ad- 
ditional COM energy Ecm of a tightly bound molecular 
state is therefore: 
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) is the single-molecule eigenstate of Hj^ 



again treated as the zero-order wave- function, and E^j^j, 

Efjj^j are the harmonic and anharmonic contributions re- 
spectively to the COM molecular energy . The reduced 
anharmonic correction compared to the free atomic case 
is due to the reduced spatial width of the wave-function, 
caused by the increased molecular mass compared to an 
atom. 
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C. Tightly bound resonance threshold 

This allows us to make a relatively simple calculation. 
A threshold condition in the tight-binding limit is ob- 
tained from equating the total ground state transverse 
energy of two atoms with the total molecular energy in an 
COM excited transverse state. It will be convenient for 
later calculations to define a dimensionless bound state 
energy relative to unbound atoms in a transverse waveg- 
uide as: 
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where rj — ujx/ujy is the anisotropy of the two transverse 
binding frequencies, and Ei, = E'^ in the strong binding 
limit of interest here. The correction term of — i (1 + 77) 
is required to take account of the difference in the trans- 
verse confinement energies between the atoms and the 
molecular state, which does not occur in free space. 

After including the anharmonic corrections and exci- 
tation energies from Eq p.Sp and Eq (|3.14p . one obtains 
a resonance condition for deeply bound molecular ACIR 
in an intuitive form as: 
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On transforming this to dimensionless form, we obtain: 
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Clearly there are multiple resonances as N^ and Ny 
is varied, thus altering the COM quantum numbers of 
the excited transverse molecular states. The position of 
these resonances is largely determined by the quantum 
numbers, together with anharmonic shifts. In the next 
section, we show that wave function symmetries mean 
that the even order resonances are directly coupled by 
the strong quartic anharmonicities ax^Uy. Odd order 



Figure 2. (Color online) Diagram of the confinement in- 
duced energy levels in an anisotropic trap or waveguide, with 
anisotropy 77. The solid line is an internally excited two- 
particle ground state. The dashed lines correspond to a COM 
excitations of the two-particle ground state. Confinement in- 
duced resonance thresholds occur when these excited molec- 
ular states become resonant with two incoming free particles 
with zero momentum. This is indicated by a star in the in- 
ternally excited CIR case, which is the only level coupled in 
the case of a harmonic trap. For anharmonic traps, these 
additional levels become coupled to input atomic states with 
zero transverse excitation, and acquire additional anharmonic 
energy shifts. 



resonances are coupled through the relatively weaker cu- 
bic anharmonic terms due to the v-Wx terms, which are 
physically caused by the slow variations in the Gaussian 
envelope function of the trapping lasers in these experi- 
ments. 

The consequences are seen in FiglU which shows the 
first two strongly coupled even-order COM resonances, as 
compared with the internally excited resonance position. 
The traditional harmonic confinement induced resonance 
(HCIR) state in a one dimensional harmonic waveguide is 
identified with the solid line in Figl^l This curve is simply 
the first internally excited state of the two-body ground 
state |l^, [13 ■ Ons expects a CIR to occur whenever 
the solid curve crosses the lowest horizontal line, thus 
permitting a resonance to occur with incoming atoms 
near zero energy. 

However, we see that there are also two further possi- 
bilities, spaced both above and below the internally ex- 
cited resonance. These are the lowest lying even order 
ACIR resonances, which we expect to be the dominant 
excitations in the case of anharmonic waveguides. 

While this calculation is approximately correct, and 
gives an excellent intuitive picture, we will carry out a 
quantitative calculation in the following sections with the 
inclusion of the internal degrees of freedom as well. 



IV. 



ANHARMONIC CIR IN A QUASI-ID 
WAVEGUIDE 



While Eq (|3.17p is valid for a relatively deeply bound 
molecule, it neglects both anharmonic and waveguide 
corrections to the internal molecular energy Eb. The 
question of which transverse molecular quantum num- 
bers Nx,Ny are accessible in terms of selection rules also 
needs to be addressed more carefully. In this section, 
we therefore treat the general case of a quasi- ID confin- 
ing waveguide with an asymmetric confining potential. 
This case can be continuously changed to the limit of the 
2D trap, which is treated in more detail in the following 
section. Due to the small anharmonic parameter a, the 
anharmonic term Ha will be treated as a perturbation to 
Hh . 



The dimensionless ground state molecular energy e is 
given by an implicit equation [19|: 



0,30 \/t^ 

where the RHS is defined by the definite integral 



(4.3) 
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In the strong binding limit, the limiting behavior of 
this integral is: 



lim 
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A. Harmonic CIR 

Firstly, we consider the anisotropic waveguide without 
the anharmonic perturbation Ha- While this has been 
treated previously [1^, we revisit it here as a first step to 
solving the anharmonic case. The origin of the confine- 
ment induced resonance is due to two scattering atoms 
forming a virtual molecule via their S-wave interaction. 
The energy of the resulting two atom quasi-bound state 
in a waveguide is written as: 



E' 



(0) _ 

M — 



EcM 



Eb 
1 



fU^'T. 



N„ 



tVjJy + Eb (4.1) 



Here Ecm is the energy of the COM excitation, {Nx,Ny) 
are the quantum numbers of the COM, and Eb is the 
binding energy of the two atoms. This becomes resonant 
with two incoming atoms near zero momentum when 
Ej^/ — Escatt, where Egcatt is the incoming free-particle 
energy. This is trivially given in the zero-momentum, 
harmonic case by: 



fUjJx + hUJy 



(4.2) 



By solving the eigenproblem of the relative Hamilto- 
nian of the two atoms HJ^'^^ , we can obtain the relation 
between the binding energy Eb and the 3D s-wave scat- 
tering length a^D . This is known from previous work [l9| , 
by solving for the dimensionless energy e of the molecular 
ground state, where e = Eb/huiy — (77 -I- 1) /2, as given in 
Eq ([XT5|) . We note that, just as with Eq ((XT5)l in the 
previous section, the dimensionless energy e is defined so 
that it includes the change in transverse confinement en- 
ergies. With this definition, Eb reduces to the free-space 
binding energy in the limit of weak confinement or strong 
binding. 



This is leads to the free-space three-dimensional bind- 
ing energy result, Eq (|3.12p . as one expects in this limit. 
Combining Eq (|3.15p and Eq (|4.ip . the threshold condi- 
tion for a harmonic trap can be summarized compactly 
in one equation as 



e + Nxi^ + Ny = Q, 



(4.6) 



so that: 



— = -^-^1 (- [N.V + Ny],Q) . (4.7) 

Further resonances are anticipated if we consider rela- 
tive atomic motion. Such an internally excited molecular 
state is described by a completely different integral equa- 
tion. For the first excited state of the internal motion, 
one obtains: 






(4.8) 



where the RHS is now defined by the definite integral 



J'e (Ce, 0) 
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The CIR threshold condition including both internal 
and COM excitations is therefore 



so that: 
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asD 



ee + N^V + Ny^O, 



-.T, (- [N^i^ + Ny\ , 0) 



(4.10) 



(4.11) 



The above equations 14. 6l and l4. 101 give threshold condi- 
tions in the limit of small anharmonicity, for coupling to 
either the internal ground or excited state respectively. 



with center-of-mass quantum numbers included in the fi- 
nal resonant state. 

As such, they give an elegant picture of the possi- 
ble resonances, including both internal and COM exci- 
tations. However, if the anharmonic term Ha is not in- 
cluded, the two incoming atoms in a transverse ground 
state cannot couple to the transverse excited molecular 
states during the collision. Hence, for harmonic confine- 
ment there is only one observable CIR resonance no mat- 
ter how anisotropic the transverse confinement [l9|. This 
is described by the last equation above, Eq (|4.11l) . on set- 
ting N^ = Ny = 0. 

In reality, anharmonic terms do occur. These lead both 
to couplings that allow COM excitations, and to energy 
shifts which alter the resonance locations. In the follow- 
ing analysis, we will assume that there are only COM 
excitations, and we will apply perturbation theory to the 
bound state COM energies predicted by Eq (|4.6p . 



B. Anisotropic, anharmonic CIR 

If the anharmonic perturbation Ha is now introduced, 
the COM motion will be mixed with the relative mo- 
tion by the anharmonicity of the confining trap. Then 
the transversely excited COM molecular states can cou- 
ple to the scattering state of the two incoming atoms in 
the transverse ground state. However, both the atomic 
and molecular states now have energy levels shifted by 
anharmonic corrections. This means that the fundamen- 
tal resonance equation is modified from Eq (|4.6p for the 
harmonic trap case. It is now: 



+ N^-q + Ny + ea {N^,Ny) = , 



(4.12) 



where ea{Nx,Ny) is the anharmonic correction to the 
relative energy levels for a COM excitation with quantum 
numbers {Nx , Ny ) . This has been treated already in the 
deeply bound limit in Eq (|3.17l) . We now treat this in 
the general case. 

From Eq (j3.8p , one must include the input anharmonic- 
ity in the atomic levels, which for the case of two atoms 
in an initial atomic transverse ground state is given by 
Eq (13. 9p . Next, we consider the effects of anharmonicity 
on the bound or molecular energies. 

With a quartic anharmonicity which is symmetric 
around the origin, there are constraints on the types of 
coupling that can occur to the COM motion. In partic- 
ular, with a symmetric input state having N^ ~ Ny ~ 0, 
one must have a symmetric resonance state. This implies 
that only even COM quantum numbers are strongly cou- 
pled. There is also a weak coupling to odd COM quantum 
numbers, caused by cubic anharmonic parameters, which 
we treat in the next section. The lowest of the strong non- 
linear resonances occur when {Nx,Ny) = (2,0) or (0,2). 
Consequently the resonance splits if the transverse con- 
finement is anisotropic. 

We give a detailed calculation of the effects of anhar- 
monic confinement on the bound-state energies of these 



resonances in the Appendix. Using these results, we ar- 
rive at the resonance condition for the N^ = 2 state , 
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In like manner, the resonance condition for the Ny 
state is, 
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These results can be compared with Eq p.l7p . which 
were obtained in the previous section, dealing with the 
case of a deeply bound molecular state. On dropping 
terms scaling with 1/e, the two conditions agree in the 
limit of strong molecular binding, with e — >■ —00. 

For experiments using optical lattices with equal wave- 
lengths in each direction, one finds that: 



axV 
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Then the resonance conditions can be reduced to a sim- 
pler form. For {Nx,Ny) — (2,0), one obtains 
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Solving these equations requires the use of a nonlinear 
equation solving numerical algorithm, like the Newton- 
Raphson method which we use in FigE] to obtain theo- 
retical results for comparison with experiment. 

Now we are ready to understand the dominant ef- 
fects observed in the recent Innsbruck experiment [23|. 
The atoms are prepared in a ID tube geometry us- 
ing two orthogonal pairs of counter-propagating laser 
fields to create a trapping environment. In the case of 
transverse isotropic confinement, due to the small colli- 
sion energy and large energy interval between the trans- 
verse energy levels, only the first excited molecular states 
{Nx =2,Ny = 0) and {N^ = 0,Ny = 2) of the COM mo- 
tion which are degenerate can couple to the scattering 
state of two incoming atoms. In this isotropic situation, 
the results of ACIR are similar to those of HCIR and 
also agree well with the experiment, as long as the an- 
harmonicity of the confinement is not too large. Due 
to the broad resonances and fitting techniques used, the 
isotropic experimental data has too little precision to ac- 
curately distinguish between ACIR and CIR resonances. 
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Figure 3. (Color online) The predictions of ACIR resonances 
using first-order anharmonic perturbation theory comparing 
with the experimental data [23l |. The black solid curve is the 
result of a generalized Olshanii's HCIR model [13|. Here the 
anisotropy rj is changed by increasing the frequency oj^ with 
a fixed ujy = 2n x f3.2(2)kHz (or ui in[23|). which resulted in 
an anharmonic parameter of a = —0.133 . 



By increasing the transverse anisotropy 77 = 
LOx/ijJy , the molecular states {Nx = 2, Ny = 0) and 
{Nx = 0, Ny = 2) of the COM motion are no longer de- 
generate. A splitting of the ACIR is expected, which has 
been observed in the experiment. By contrast, no split- 
ting of the linear HCIR resonance is predicted. Thus, 
the double resonance observed experimentally is a clear 
signature of the ACIR resonance. It is also a sign that if 
there is any internally excited HCIR effect in this exper- 
iment, it is relatively suppressed compared to the ACIR 
effect which involves COM excitations. 

The theoretical predictions of the CIRs compared with 
the experimental results [23| are presented in FiglSl As 
we can see, the predictions of the COM-relative coupling 
theory are consistent with the experimental data. 

In these experiments, the resonances were identified as 
occurring at the point of maximum molecular loss. This 
leads to an offset between the apparent and true reso- 
nance due to finite resonance widths. To compensate for 
this in the original experimental publication, the data 
was fitted to the isotropic HCIR prediction by adding a 
small constant offset to make it equal to the theoretical 
results at 77 = 1. We follow a similar fitting procedure 
here, for consistency, but we fit the offset to the isotropic 
ACIR prediction instead. Note that the anharmonic pa- 
rameter is different to that used in TableU because Table 
|T] gives values representative of a range of experiments, 
while in each of the figures we use the data from the 
relevant experiment. 

At the point 77 = 1, there is transversely isotropic con- 
finement, and the ACIR prediction is in accord with the 
Olshanii CIR, apart from anharmonic corrections. How- 



ever, as the transverse confinement becomes more and 
more anisotropic, the internally excited HCIR resonance 
persists as a single resonance except for a small frequency 
shift [19]. Both ACIR theory and experiment show a 
clear splitting of the original resonance with increasing 
trap asymmetry. There is excellent quantitative agree- 
ment in the amount of splitting. In these experiments, 
there is strong evidence for the COM excited ACIR res- 
onances. 



C. Multiple resonance ACIR at large anisotropy 

As the transverse anisotropy 77 becomes even larger, 
the energy spacing in the y direction decreases, and more 
transversely excited molecular states can readily couple 
to the initial scattering state. Consequently, more addi- 
tional COM resonances can occur. By contrast, there is 
still only one internal HCIR resonance predicted no mat- 
ter how large the transverse anisotropy. However, several 
multiple resonances are observed experimentally at large 
anisotropy. In order to understand these, we now con- 
sider the other transverse states. 

As in the previous section, given any coupling param- 
eter (asij), the unperturbed binding energy e is deter- 
mined by Eq (|4.3p . which as an implicit equation in- 
volving the integral Fi (e,0). We then use perturbation 
theory to calculate the dimensionless anharmonic energy 
shift Cq. Hence, we obtain the ACIR resonance position 
of osD for arbitrary odd and even COM quantum num- 
bers from the solutions to the overall resonance equation, 
Eq Km . 

For brevity, we will refer to an arbitrary molecular 
COM resonance as simply {Nx,Ny). The general form 
of the resulting anharmonic shifts is derived in the Ap- 
pendix, for arbitrary quartic anharmonic parameters. We 
include both odd and even order resonances, because, as 
remarked earlier, there are both cubic and quartic anhar- 
monic couplings, which leads to the possibility of both 
even and odd ACIR resonances. However, for simplicity 
we do not include the relatively small cubic energy shifts. 

For comparison to current experiments, we are inter- 
ested in the case of equal optical trapping wavelengths, 
which means that axi] = ay = a. From the Appendix, 
the general form of the resulting anharmonic shifts for 
resonances corresponding to the {Nx,Ny) COM state is 
as follows, 



eaiNx,Ny)/ 



3iN^ + Ny + N^ + Nx)-9 
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This allows the positions of ACIR denoted by 0313 to be 
calculated. In the experimental reports of multiple reso- 
nances at large anisotropy (Fig. 4 [23j), the loss rates are 
plotted as a function of magnetic field, rather than a^D ■ 
We therefore make use of the relation between the 3D 
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scattering length asnand -B-field for the relevant ^^''^Cs 
Feshbach resonance [3ll|, which is, 



and 
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(4.20) 



Hence, the predicted magnetic field at resonance can 
be calculated. In order to compare our theory to these 
experiments, the regime of anisotropy ry considered is 
[1.4,2.3] . We also note that in these experiments, the 
trapping frequency LOy is varied, so the effective an- 
harmonicity parameter a therefore changes at different 
anisotropy. 

Then main results are summarized in FigUl which 
compares theory to experiment. The experimental reso- 
nance points are obtained from the raw data as the start 
of the resonance edges, which is appropriate for this type 
of resonance experiment. No fitting parameters or shifts 
are employed. Since error bars were not given in the 
experimental plot, we are unable to estimate these. 

Most observed resonances can be easily identified, 
which are coded in the same color as the correspond- 
ing theoretical predictions. There are 25 identified reso- 
nances, all of which are in excellent agreement with the- 
oretical calculations. However, there are 5 smaller reso- 
nances not identified, which are indicated by the empty 
triangles. 

Apart from experimental issues, possible explanations 
for these unidentified peaks are: 

• Higher-order many-body ACIR resonances. Com- 
bination 4-body cluster resonances could occur at 
intermediate points between the identified 2-body 
resonances. For example, the three unidentified res- 
onances between (0, 1) and (0, 2) could be caused 
by the simultaneous excitation of (0, 1) and (0, 2) in 
a 4-body collision. Similarly, the three resonances 
we have identified as (1,1) resonances could also 
be caused by 4-body excitation of (0, 3) and (0, 4) 
ACIR resonances. 

• Anharmonically shifted internal HCIR resonances. 
We have not calculated these, as the anharmonic 
shifts of these resonances are outside the scope of 
this paper. However, this mechanism provides a 
possible explanation for the two unidentified reso- 
nances between (0,4) and (2,0). 



ANHARMONIC CIR IN A QUASI-2D 

SYSTEM 



For a quasi-2D system, atoms are tightly confined in 
the axial direction, and can freely move in two trans- 
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Figure 4. (Color online) The resonant magnetic field B vs 
anisotropy r\ at ACIR resonances (solid and dotted lines iden- 
tified alongside the figure) from first-order anharmonic per- 
turbation theory comparing with the Innsbruck experiment 
data (round circles). Data is extracted from multiple reso- 
nance scans [23 | at absorption edges. Here the anisotropy 
r\ is changed by decreasing the frequency lOy with a fixed 
oja; = 27r X 16.6(2)kHz which results in an 7;-dependent an- 
harmonic parameter a = — 0.106?7 . 



verse directions. In the following section, we will derive 
the anharmonic CIR resonance properties in this case as 
well. As we will show, these results can be also obtained 
from results in the previous section by taking one of the 
confinement frequencies to zero. However, the direct cal- 
culation given here is an important check on the consis- 
tency of our approach. In the case of an optical lattice, 
the trapping potential is of form: 



[/«^* (r) = K cos^ 



27rz 



(5.1) 



mij?z^ ( 1 -f 



az 
~d? 



Here, Vz is the potential well-depth with one standing- 
wave trapping laser at optical wavelength A . This leads 
to the trap frequency uj in the z direction. We note that 
cubic anharmonic terms are also possible due to focusing 
effects in this case also. 

As before, we use a scale length of the reduced oscil- 
lator length d — ^J2h/rriuj , and define a single atom 
oscillator length a±^ = y^h/muj, and a small anharmonic 
parameter a, so that for an optical lattice. 






(5.2) 



The Hamiltonian of two atoms in a quasi-2D system is 
then, 



H^Hi+H2 + U{vi~ rs) , 



(5.3) 
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where, 



The anharmonic energy Ej^J 



(5.4) 
As previously, U (jci — r2) is the interatomic interaction, 
and we set 

R=^(ri+r2) 
r = ri r2 . 

Then Ea. (|5.3p becomes, 

H = Hh + Ha 
Where, 



H,, = - 



2M 



R 



-Mu^Z' 



(5.5) 



(5.6) 
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(1) 



Including the first order modification _B^^ 
at the following integral equation for e. 



we arrive 



N - 



hid 



hjj + E\j (e) — hoj + —ahoj 

(5.12) 

By solving this equation, the binding energy e*-^-* at the 
resonance is obtained. Then substituting e'^-'^^ into the 



'■3D 



at 



Eq. (|5.10p , we will obtain the 3D scattering length 

the resonance. Following the general procedure outlined 

in the Appendix, we find that a resonance occurs at: 
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Hence, we obtain the equation that e should satisfy for 
states with N = 2 , 
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Here r = (x,?/, z) , R = {X,Y,Z) , and M = 2to , 

/i = ?7l/2 . 

Due to the small anharmonic parameter a , the anhar- 
monic term Ha can be treated as a perturbation to Hh . 
Firstly, consider the case without the small perturbation 
Ha ■ An anharmonic confinement-induced resonance is 
expected when the energy of the virtual molecule is de- 
generate with the energy of two incoming atoms from the 
axial ground state, 
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Where Ecm = (-^ + 5) ^^; Egcatt = So;, and the binding 
energy Ejy is determined by. 



d 
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and the integral expression required here is given by, 

exp (^ei) 1 



^2(e,0) 



dt 



tVT 



t3/2 



(5.11) 



We define Ei, = {e + 1/2) fwj , d = y^h/jiuj . The lowest 
resonance occurs as TV = 2; however, this is not coupled 
to the incoming states unless there is an anharmonic term 
in the potential. 

If the anharmonic term Ha is included, we need to 
consider how the energy of the input atomic states and 
virtual molecule Em is affected by this perturbation. 



This result is identical to Eq (|4.14p in the limit of 
?7 — > 0, as one might expect, since the quasi-lD trap 
becomes two-dimensional in this limit. However, it is in- 
structive that one cannot regain this limit directly from 
Eq (14.171) . which holds for optical lattices with equal 
wavelengths. The reason is very simple: in an optical lat- 
tice at low transverse confinement frequency, the trans- 
verse wave-function becomes more and more deconfined. 
This increases the relative anharmonicity, as given in Eq 
p.4[) . so that a ^ l/w. Therefore, our anharmonic per- 
turbation theory would break down for a weakly con- 
fined ID system described by Eq (|4.17p . Optical lattices 
that are only weakly confining require a full Bloch wave- 
function theory, typically requiring detailed numerical di- 
agonalization |32| . 

From Eq. (l5.14ll . we can calculate the relative binding 
energy e for the states of iV = 2, in the two-dimensional 
limit. Then, substituting into Eq (|5.10p . the correspond- 
ing resonance scattering length a^jj is obtained. Here, 
in order to calculate e numerically, we use the Newton- 
Raphson method as before, together with a numerical 
calculation of J^2 (e^O). The position of the ACIR de- 
noted by the 3D scattering length as a function of the 
anharmonic parameter ad'^ is presented in FigjSl and we 
can see that the 2D ACIR is predicted at the regime 
of asD > . This is in contrast to the internal HCIR 
resonance, which occurs for the attractive regime with 
asD < 0. 

Now we can compare these predictions with 2D reso- 
nances observed in the recent Innsbruck experiment [23j. 
When one of the lattice lasers is turned off, the system 
approaches a 2D geometry. All resonances disappear ex- 
cept one at a^D > . Given the average anharmonicity 
in this experiment of a = —0.121, we expect as shown in 
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Figure 5. (Color online) Predicted resonance positions of 2D 
CIRs, as a function of the dimensionless anhamionic param- 
eter a. The experimental data point is an average from ex- 
periments [231, plotted at the average anharmonicity value. 
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Figure 6. (Color online) Predicted resonance positions of 2D 
ACIR resonances, as a function of the transverse confinement 
parameter, a±. The experimental data points are taken from 
experiment [2a |. where circles indicate resonances identified 
by absorption edges. 



FiglS] to find a resonance at 



0-30 = 0.85aj^ 



(5.15) 



The observed resonances occur at a constant ratio be- 
tween asD and a±, such that 



^^asD = 0.84(3) 
a± 



(5.16) 



Tliis is consistent with the prediction of 2D ACIR. The 
anharmonicity a in the experiments was varied through 
a small range as the trapping frequency varied, and we 
have calculated the ratio at the average anharmonicity. 

Finally, we plot the predicted variation of a^^ with 
transverse confinement parameter a j^ , at a fixed value of 
ua corresponding to the Innsbruck experiments, in Fig|51 
This data is also in excellent quantitative agreement with 
ACIR predicted resonance positions. 

In anhamionic trap experiments, we would generically 
expect both the ACIR resonance at a^o > 0, and the in- 
ternal HCIR resonance at a^n < to be observed. The 
Innsbruck experiments [23] show clear evidence for the 
ACIR resonance, but not the internal HCIR resonance. 
The Cambridge experiments [2J| have shown evidence for 
HCIR resonances, although the magnetic field scans were 
not large enough to observe any ACIR resonances. The 
question of which is observed will depend on the size of 
the anharmonic parameter, the dynamics of the experi- 
ment, and the method of detection of the resonance. It 
appears likely that the molecular loss technique is par- 
ticularly sensitive to COM resonances. 



VI. CONCLUSION 

We have extended the theory of the harmonic 
confinement-induced resonance to include the effects of 
anharmonic confinement, or ACIR, since previous har- 
monic theories cannot explain phenomena observed in 
recent experiments. In the presence of anharmonic per- 
turbation of the confinement trap, the COM motion of 
two atoms couples to the relative motion, and additional 
resonances appear. We have calculated the energy of the 
resulting new resonances up to first order in perturba- 
tion theory. The results agree well witli experiments, 
with both even and odd order multiple resonances being 
found. These differences are not just small perturbations 
on previous HCIR predictions, and show large qualitative 
differences from the predictions for harmonic traps. 

ACIR resonances due to COM excitations are always 
present in Feshbach systems with any form of transverse 
confinement. The important issue is that for harmonic 
traps the Kohn theorem prevents these resonances from 
being coupled to incoming atoms in the two-body sector 
of the COM ground state. However, experimental optical 
traps do have relatively large anharmonic parameters, 
which allows the resonances due to COM motion to be 
coupled to incoming scattering states. We find excellent 
quantitative agreement between ACIR predictions 
and recent experimental observations of confinement 
resonances, both in one and two dimensional traps. 
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Appendix: Anharmonic Perturbation Theory 



We consider how the relative anharmonic energy of the virtual molecule Ca Eq. (|4.12[) , is calculated from the anhar- 
monic perturbation Ha, including the effects of anharmonicity on both relative and COM motion. For space reasons, 
we do not give all the overlap integrals, as the calculations are similar in all cases. Instead, we focus on the most 
important cases. 

First, note that the dimensionless change in the relative anharmonic perturbation energy is: 



1 



p(l) 9p(l) 



(A.l) 
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(1) 



where E\ is the perturbed atomic anharmonic energy given in Eq p.Sp , and E)J is the first order perturbation 
theory correction to the energy of the virtual molecule: 



j(i) _ 



EVJ = (^AT^AT^bl Ha IV'&^AT.wJ • 



(A.2) 



Here. N? 



N:^N: 



M^b) 



* 



N^N„ 



\ipi,) is the eigenstate of Hh- This is treated as the zero-order wave-function, and 



includes both relative and COM motion. Next, we introduce 
Hamiltonian, so that 



with the standard two-dimensional harmonic oscillator solution of: 

exp (-2e^,,) il^,,^ (2C.,,) 



^jv^AT ) as the unperturbed wave-function of the COM 

(A.3) 



»Wx 



,) 



Vny^d^^y2^^.y-^N,J 



(A.4) 



p(i) 



where £^i — {X/dx, Y/dy). We now turn to the task of evaluating Ej^/ , the anharmonic correction to the total molecular 
energy at a given unperturbed dimensionless energy e. 



1. The relative wave-function \tf>b) 

The specific form of \tph) can be obtained by directly solving the eigenproblem of Hf-^^^ , 



^i^ 1 1 

-^V2 + -^iu;lx^ + -fiujy + V (r) 



= E\ 



The wave-function 1-0) can be written as. 



IV') = ^0 (r) - / dr'GE (r, r') V (r') ^ (r') 

, , 2Trh?a-ir) „ , 
= Vo(r)-7^ —GEir,0) 



Where tpo (r) is the regular solution of Ea. (|A.5P and, 

n = 



—r^ (r) 



(A.5) 



(A.6) 



(A.7) 
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Here, we have use the pseudo-potential approximation. For a bound state, the solution of Eq. (jA.5l) 1-0) should vanish 
as r — ^ CXI, hence ipQ (r) = . Then the constant TZ is only a normalization coefficient. The Green's function Ge (r, r') 
is the solution of the following equation. 



^ 2 2,-'^ 2 2 



Ky - E 



Ge ir,r')^ 6 (r-r') 



2^ "^ 2^ ^ 2' 

The form of the Green's function Ge (r, 0) for a bound state can be easily calculated. 



Ge<o(r,0) 



V^-^i ("") „„„ f ^^^ + y^ 



27r3/^d^hu}, 



■ exp 



2(i2 



where we introduce the integral representation: 



/i(r) 



dt- 



explf^-lS 



e ' r]x 



e y 



10 ^t(l-e-"*)(l-e-*) 

Recalling that rj — Lu^/ujy , and E = {e + rj/2 + 1/2) ftwj,, the bound-state wave-function is, 



\ipb) =A/'exp 



2d2 



/fr) 



where A/" is the normalization coefficient. 



2. ID COM overlap terms 



(A., 



(A.9) 



(A.IO) 



(A.11) 



We can now calculate the overlap integrals which give the anharmonic contributions to the quasi-lD and quasi-2D 
ACIR resonances. In order to calculate the first order modification Ej^.^ , we focus on the {Nx,Ny) COM state. The 
relevant terms are: (*jv^JvJ^^*Ar,Ar„) , {^N^Ny\y^ I^JVxivJ , (*JV.JvJ ^^ I'^n^n^) , {^n^nJY"^ I'^N^Ny)- 

We define D^^y — y''h/2miOx^y = dx^y/2 as the transverse confinement parameter of an atom pair in the following 
calculation. Hence: 

2 



{'^N^Ny\X^\ipN^Ny) 



^Dx2^^NJ 

1 

V^Dx2N^NJ 
N. + l 



X^exp(-I) 



oo 

5 rzoN, 



■ DlV^2^^NJ 



N,. 



1 
"2 



X 



dX 



N^ + l) + l^x i^x - 1) 



N- + l) + 1^- (^- " 1) 



£>■* 



(A.12) 



(V'. 



'^^^y\^'\^^^^y)= ^Dx2^^NJ 



OO 



X^exp(-|; 



oc 

,3 ^r,N, 



■ Dly^2^-NJ N, 



^Dx2^^Nx\ 

By using the exchange symmetry of x and y , we can easily obtain, 

1\ f .. 3 



{'^N^Ny\Y'^\tpN^Ny) 



N„ 



{-^N^NylY^ l^N^Ny) ^ [Ny + lj 



Dl 



Here, we have used the formulas, 



/oo 
dtt^e-*^Hl (t) = v^2"n! 
-OO 



dtt^e^''H^^{t) = V^2''nl 



X 



Hm — 



dX 



Ny + 7;]+7;Ny {Ny - 1) 



y 



n+- 



3\ 1 , 

— H — n n 
2/ 2 ^ 



(A.13) 

(A.14) 
(A.15) 

(A.16) 
(A.17) 
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3. ID Internal motion overlap terms 



The relevant terms are now: {ipt,\ x^ \il)b) , (V'&l V* \i^b), (V'f;! ^^ IV'fe) i (V'fel V^ {''Pb)- For the bound state of the relative 
motion \ipi,) , the main contribution of the integrals, e.g., {ijji,\ x^ j'i/'b) etc., comes from the regime around the origin. 
We note that the energy parameter e is negative. The form of j'i/'b) near the origin r = can be easily obtained from 
Ea. lTOTI) . 



1 exp I 



-2€r/dy 



/TuL, 



which has been normalized. Then, 

„2 



(V'bla;^ IV'b) = — i-\/-r- / / / r^ sin^ 6* cos^ ^ exp ( ] drd9d(j) 

Trdy V 2 Jo Jo io V "" 



Xj, 



12e 



(A.18) 



(A.19) 



2'7r /'TT /'C50 



(^-61x4 1^6) = ^./-^/ / / r^sin^ycos 
5dt 



4 exp [-2^/-2er/dy) ^ 



r sin OdrdOdcf) 



40e^ 



(A.20) 



1 ^^^^ /'27r /'TT /'OO 

■nay \ z jQ Jo Jo 

= -^ 
12e 



r^ sin^ 61 sin^ (t>^^ „ '^^' ^K ^ sin edrd9d(t> 



(A.21) 



1 / /'27r /'TT /'OO 

TTflj, V 'i Jo Jo Jo 

3dt 



r^ sin'' e sin'' </.^^^^^-^^ y__^!j_Jii^2 ^j^^ Qdrdedc/) 



40e2 



(A.22) 



4. The first order modification of the energy E 



(1) 



Using the overlap integral results given above, we find that: 



(i^N^Nyi^b] HI \i)N^Ny^b) 



32 



(V'WxW^&l Hi \%l^N^Ny^b) = ^g^ ^ 



6A^x (AT, + 1) + 3 - 



6iVj,(A^y + l) + 3 



r;(2A^, + l) , 3772 



(2A^y + 1) 



10e2 
3 

l(k2 



(A.23) 



(A.24) 



Combining this with the atomic energy correction gives the overall result for an equal wavelength ID optical lattice, 
in which ax^jJx — otyUJy : 



.(!£') 



32 



6{N^ +Ny + Nl + Nx) - 18 - 



77 (2A^, + 1) + (2A^y + 1) , 3(7?2 + l) 



e 10e2 

Following similar procedures, from Eg IA.241 the anharmonic correction in the corresponding 2D case leads to 

(2Ar + l) 3 



(A.25) 



,(2£>) ^ ^ 
32 



6A^(iV + l) -9- 

r 



10e2 



(A.26) 
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